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m ■ ABSTRACT 

^ ■ We present models of photometric evolution of galaxies in which the effects of a dusty 

OO . interstellar medium have been included with particular care. A chemical evolution code 

^ i follows the star formation rate, the gas fraction and the metallicity, basic ingredients 

'- ^ ' for the stellar population synthesis. The latter is performed with a grid of integrated 

■ ■ spectra of simple stellar populations (SSP) of different ages and metallicities, in which 

the effects of dusty envelopes around asymptotic giant branch (AGB) stars are included. 
^ I The residual fraction of gas in the galaxy is divided into two phases: the star forming 

I molecular clouds and the diffuse medium. The relative amount is a model parameter. 

^ ■ The molecular gas is sub-divided into clouds of given mass and radius: it is supposed that 

^ . each SSP is born within the cloud and progressively escapes it. The emitted spectrum 

^ I of the star forming molecular clouds is computed with a radiative transfer code. The 

■ - - ' diffuse dust emission (cirrus) is derived by describing the galaxy as an axially symmetric 

system, in which the local dust emissivity is consistently calculated as a function of the 
local field intensity due to the stellar component. Effects of very small grains, subject 
to temperature fluctuations, as well as polycyclic aromatic hydrocarbons (PAH) are 
included. 

The model is compared and calibrated with available data of normal and starburst 
galaxies in the local universe, in particular new broad-band and spectroscopic ISO ob- 
servations. It will be a powerful tool to investigate the star formation, the initial mass 
function (IMF), supernovae rate (SNR) in nearby starbursts and normal galaxies, as 
well as to predict the evolution of luminosity functions of different types of galaxies at 
wavelengths covering four decades. 

Subject headings: dust, extinction — galaxies: ISM — galaxies: spiral — 
galaxies: starburst — infrared: galaxies — radiative transfer 
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1. Introduction 

Observations performed in tlie last decade or so, 
particularly in the IR regime, have clearly demon- 
strated that dust is one of the most important com- 
ponents of the interstellar medium (ISM), containing 
a large fraction of heavy elements ejected from stars. 
This paper is mainly concerned with the influence of 
dust grains on the transfer of radiation emitted by 
stellar systems. Basically dust absorbs and scatters 
photons, mostly at wavelengths ;$ 1 /xm and returns 
to the radiation field the subtracted energy in the 
form of IR photons. The resulting spectral energy 
distribution (SED) is often substantially changed and 
in many relevant cases radically modified. 

Not surprisingly, dust reprocessing of the optical- 
UV photons emitted by stars into the infrared regime 
turns out to be particularly severe in galaxies under- 
going massive episodes of star formation, which pref- 
erentially occur within dense molecular clouds. In- 
deed dust could affect galaxy evolution because it 
modifies the physical and chemical conditions of the 
ISM. In particular star formation is at least favored 
by the presence of dust, which shields dense clouds 
from stellar UV radiation and keeps them to temper- 
atures low enough to allow the onset of gravitational 
instability. 

The SED of normal star-forming disk-like galax- 
ies is also affected by dust, mainly associated with the 
diffuse ISM and with the envelopes of evolved asymp- 
totic giant branch (AGB) stars. As for early type 
systems, it has been often suggested that the first 
episode of massive star-formation could be essentially 
hidden to optical searc;hes due to dust reproc;c;ssing, 
since their ISM might have been metal enriched on 
a very short timescale, the lifetime of the first gen- 
eration of massive stars (Franceschini ct al. 1994 and 
references therein; Cimatti et al. 1997). Even local 
quiescent elliptical galaxies exhibit infrared emission 
(e.g. Jura 1986; Bally & Thronson 1989; Knapp et al. 
1989; Roberts et al. 1991) and visual dust obscuration 
(e.g. Veron-Cetty fc Veron 1988; Kim 1989; Goud- 
frooij et al. 1994), arising in part from dust grains 
continuosly formed in dusty outflows from AGB stars 
(Tsai & Mathews 1996; Bressan, Granato, & Silva 
1998), as well as molecular line omission (e.g. Gordon 
1990; Roberts et al. 1991; Lees et al. 1991). Over- 
all, a direct comparison of the luminosity functions of 
galaxies in the optical and in the far-IR shows that, 
locally, about 30% of starlight is dust-reprocessed. 



At least three different dusty environments must 
be taken into account in order to properly understand 
the UV to sub-mm properties of galaxies: (i) dust in 
interstellar H I clouds heated by the general interstel- 
lar radiation field (ISRF) of the galaxy (the 'cirrus' 
component), (ii) dust associated with star forming 
molecular clouds and H II regions and (iii) circum- 
stellar dust shells produced by the windy final stages 
of stellar evolution. These environments have differ- 
ent importance in different galactic systems at various 
evolutionary stages. 

Despite this, in many papers dealing with the 
spectrophotometric evolution of galaxies the radiative 
processes occurring in a dusty interstellar medium 
were originally neglected (e.g. Bruzual & Chariot 1993 
and Leitherer & Heckman 1995, the former the most 
adopted one to interpret data on different galaxy 
types, the latter appositely built for starburst galax- 
ies). Indeed their codes are nowadays commonly 
used together with some prescription to approximate 
the effects of the ISM. In other cases (e.g. Guider- 
doni & Rocca-Volmerange 1987; Mazzei, De Zotti, 
& Xu 1994; Langon & Rocca-Volmerange 1996; Fioc 
& Rocca-Volmerange 1997), the effects of dust are 
included only partially and/or with substantial sim- 
plifications. For instance often extinction is consid- 
ered but thermal reradiation is not, or scattering is 
neglected, or not all the relevant dusty environments 
are considered, or optically thin emission is assumed, 
or a unrealistic geometry is adopted. Those works 
facing a complete computation of the radiative trans- 
fer through a dusty medium, do not include it in the 
more general framework of galaxy evolution, but are 
instead interested either in the interpretation of sin- 
gle objects (e.g. Kriigel & Siebenmorgen 1994) or in 
understanding the effects of dust geometry on the 
extinction of the emerging spectrum (e.g. Bruzual, 
Magris, & Calvet 1988; Witt, Thronson, & Capuano 
1992; Cimatti et al. 1997; Wise & Silva 1996; Bianchi, 
Ferrara, & Giovanardi 1996; Gordon, Calzetti, & Witt 
1997). 

The inclusion of dust effects into spectrophoto- 
metric codes leads to many difficulties. While the 
integrated photometric properties of an hypotheti- 
cal dust free galaxy are geometry independent, dust 
introduces a strong dependence on the distribution 
of both stars and ISM. Moreover the optical prop- 
erties of dust grains and their dependence on en- 
vironmental conditions have not yet been fully ex- 
plored and understood. The ensuing uncertainty may 



2 



only be parametrized to some extent. The effects of 
scattering, absorption and emission of grains, which 
may not be in radiative equihbrium with the radi- 
ation field, complicate the solution of the integro- 
differential transfer equation in a complex geometry. 

Nonetheless a more realistic and complete compu- 
tation of the spectrophotometric galaxy models in- 
cluding dust is required, since neglecting the com- 
plexity of the dust effects can lead to erroneous esti- 
mates of many interesting quantities such as the star 
formation rate (SFR). Also, the estimate of the age 
of a galaxy through the fit of its broadband SED is 
hampered by the degeneracy between the colors of an 
old galaxy and those of an extinguished young galaxy 
(e.g. Cimatti et al. 1997), further complicating the 
well known age-metallicity degeneracy. 

We have therefore developed chemo-spectrophotometric 
self-consistent galactic models including dust. Three 
different dusty environments are considered: envelopes 
of AGB stars, molecular star-forming clouds and dif- 
fuse ISM; the dust model is comprehensive of nor- 
mal big grains, small thermally fluctuating grains and 
polycyclic aromatic hydrocarbons (PAH) molecules; 
for each grain family, the appropriate computation of 
absorption and reemission is performed. The models 
are suited to simulate galaxies of any Hubble type, 
in different evolutionary stages, since we include the 
possibility of different geometrical distributions both 
for stars and gas. The complete radiative transfer 
equation is numerically solved whenever necessary. 

The many parameters and uncertainties introduced 
by the presence of dust can be constrained only by 
means of the multiband approach pursued in this 
paper, where the spectrophotometric properties of 
galaxies are consistently reproduced from the UV to 
the sub-mm, i.e. taking into account extinction at 
short wavelengths and the consequent thermal emis- 
sion in the IR regime. The combination of observa- 
tions from HST, Keck, ISO and ground-based optical, 
IR and sub-mm telescopes already provided several 
objects at significant redshift with spectral informa- 
tion on this large A range. Their number will increase 
after the completion of ISO surveys, and will burst 
when SIRFT, FIRST, Planck Surveyor and NGST 
will operate. 

In this paper we present the model and we use it to 
study local galaxies. In a forthcoming paper we will 
exploit it to study the cosmic evolution of different 
types of galaxies. 



2. Model description 

We divide the problem of estimating the SED of a 
galaxy at age tc, including dust effects, in two steps: 

(1) the history of the star formation rate ^'(t), of the 
initial mass function (IMF), of the metallicity Z{t) 
and of the residual gas fraction g{t) is determined and 
then (2) the integrated SED of the galaxy is predicted 
taking into account all the stars and the gas present 
at to. When the effects of dust are neglected, step 

(2) simply involves a sum of all the spectra of stars. 
In our case we have instead to introduce a specific 
geometry for gas and stars, and then to compute the 
radiative transfer of the radiation emitted by the stars 
and the dust. 

The purpose of the present work is mainly to de- 
velop a general procedure for step (2). Step (1) in- 
cludes a number of possible choices which can be eas- 
ily interfaced to our code in order to describe the evo- 
lution of the SED. 

2.1. Chemical evolution 

The chemical evolution is a preliminary process in 
our code. We siimmarize here the main aspects of the 
code, which basically follows the guidelines given by 
Tantalo et al. (1996). 

The code describes one-zone open models includ- 
ing the infall of primordial gas {M{t) oc exp(— t/tj„/)), 
in order to simulate the collapse phase of galaxy for- 
mation and, when required by the astrophysical situ- 
ation under study, galactic winds. 

The adopted SFR is a Schmidt-type law, i.e. pro- 
portional to some power (between 1 and 2) of the 
available gas mass: ^{t) = vMg{t)'^. For starburst 
galaxies we add to the general smooth SFR one or 
more bursts of star formation. As for the IMF we 
used the usual Salpeter law: $(M) ex M^^, with 
X — 2.35. The upper limit of the mass range is fixed 
to 100 Mq. 

This kind of chemical evolution models has been 
successfully tested against nearby spheroidal galax- 
ies, and also well reproduces the properties of differ- 
ent types of galaxies with appropriate choices of the 
parameters (e.g. a quiet spiral-type evolution is well 
mimicked with a low value of v and/or a long infall 
time scale, see Matteucci 1996 for a thorough review). 
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2.2. Synthesis of starlight spectrum (includ- 
ing dusty envelopes of AGB stars) 

The library of isochroncs for the simple stellar pop- 
ulations (SSP), the building-blocks of galaxy mod- 
els, is based on the Padua stellar models (Bertelli 
et al. 1994), with a major difference, consisting in 
the computation of the effects of dusty envelopes 
around AGB stars. The SSPs span a wide range 
in age, from 1 Myr to 20 Gyr, and in metallicity, 
Z = 0.004,0.008,0.02,0.05,0.1, to realistically repro- 
duce the mix of age and composition of the stellar 
content of galaxies. The inclusion of a wide range 
of metallicity is imposed by the aim of comparing 
model estimates to the excellent data available even 
for young high redshift galaxies. 

The spectral synthesis technique, for the starlight 
alone, consists in summing up the spectra of each stel- 
lar generation, provided by the SSP of the appropriate 
age and metallicity, weighted by the SFR at time of 
the stars birth (e.g. Bressan, Chiosi, & Fagotto 1994). 

Effects of dust in the envelopes of evolved stars 
are usually neglected in the synthesis of a compos- 
ite population. This is not justified for intermediate 
age population clusters, whose brightest tracers are 
AGB stars, and/or in a wide-wavelength synthesis 
approach. To overcome this limitation, we computed 
new isochrones and SSPs in which, along the AGB, 
a suitable dusty envelope is assumed to surround the 
star. For this envelope the radiative transfer is solved 
by means of the code described by Granato & Danese 
(1994). The envelope parameters, in particular its 
optical thickness, which is straightforwardly linked 
to the mass-loss rate and to the expansion velocity, 
are derived as a function of basic stellar parameters 
(mass, luminosity, radius, and metallicity) combining 
hydrodynamic model results and empirical relations. 
A detailed description of the adopted procedure can 
be found in Bressan et al. (1998). 

2.3. Geometry 

As already pointed out, the dust effects on SEDs 
depend on the relative distribution of stars and dust. 
The assumed geometry is sketched in Fig. 0. The 
galaxy is approximated as a system having azimuthal 
symmetry as well as planar symmetry with respect to 
the equatorial plane. We take into account three com- 
ponents: i) star forming molecular clouds complexes, 
hereafter MCs, comprising dusty gas in a dense phase, 
H II regions and very young stars embedded in it 



(young stellar objects YSO); ii) stars already escaped 
from these dense clouds (henceforth /ree stars); iii) 
diffuse gas (cirrus). 

We work in spherical coordinates {r,9,(p). The 
densities of the three components pmc, P*, and pc 
respectively, depend on r and 6 through analytical 
laws. In the following, unless otherwise explicitly 
stated, pmc and p^ have identical spatial dependence. 
In order to describe disk-like galaxies, we use a dou- 
ble exponential of the distance from the polar axis 
i? = r sin 6 and from the equatorial plane z — r cos 6: 

p^ Po e-Kp{-R/Rd) exp(-|z|/z<j) . (1) 

In the code the scale lengths Rd and Zd can be inde- 
pendently set for the three components, but in the 
models presented in this paper (Sec. ^ we simply 
adopt identical values for them. In other words we 
use two adjustable parameters when describing disk- 
like galaxies: R*^ = R'^" = i?^ and z*^ = z™'^ = z%. 
Observations within our Galaxy suggest ^ 0.35 
kpc. As a reference the e-folding scale length is re- 
lated to the absolute magnitude by 

logio(i?d/kpc) ~ -0.2Ms - 3.45 , (2) 

(Im et al. 1995). 

In the case of spheroidal systems we adopt for both 
stars and dust spherical symmetric distributions with 
King profile: 

P = Po{l + {r/r,f)-\ (3) 

extended up to the tidal radius rt . This truncation ra- 
dius is required in King models since M(r), the mass 
contained within r, diverges as r ^ cxd. Our results 
are not very sensitive to its precise choice, thus we 
simply adopted the standard value log(rt/rc) = 2.2. 
As for the stellar component we simply set 7 = 3/2. It 
has been suggested that the core radius Tc correlates 
with the luminosity (Bingelli, Sandage, & Tarenghi 
1984): 

( rc\ f -0.3 [Mb + 22.45) if Mb < -20 , 
\kpc ; 1 -0.1 {Mb + 27.34) if Mb > -20 . 

(4) 

The distribution of dust in spheroidal systems is 
poorly known, but it has been suggested that pdust oc 
/9^t„, with n ~ 1/2-M/3 (Froehhch 1982; Witt et al. 
1992; Wise & Silva 1996), i.e. 7 ~ 0.5-^-0.75. In other 
words the dust distribution seems to be less concen- 
trated than that of stars. In particular an additional 
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cool component located in the outer parts has been 
invoked to explain IRAS observations (Tsai & Math- 
ews 1996). 

The volume emissivity (ergs cm~^ s^^ stcr^^ A~^) 
of the galaxy at each point is the sum of three terms, 
arising respectively from the three components listed 
above: 

Ja^jT+jI+Ja- (5) 

In the following sections we will describe how these 
quantities, as well as the specific flux measured by 
an external observer, are computed. A preliminary 
stride is to define the optical properties of dust. 

2.4. The dust model 

Once the geometrical arrangement is specified, the 
effects of dust on radiative transfer depend on physi- 
cal and chemical properties of grains, which affect the 
way they absorb and emit photons. It is expected and 
observationally well established that these properties 
are functions of the particular environment in which 
grains happen to live. Several populations may be 
distinguished (see Dorschner & Henning 1995 for a 
review ): (i) stellar outflow dust (further subdivided 
into the two chemical subgroups of carbon rich and 
oxygen rich dust); (ii) dust in the diffuse ISM; (iii) 
dust in molecular clouds and dust around YSOs. As 
for stellar outflow dust, we refer to the treatment of 
Bressan et al. (1998) mentioned above. Here we dis- 
cuss our choices for the other two dusty components 
of the model galaxy. 

Although a lot of efforts have been devoted to de- 
rive a so called standard model for dust in the dif- 
fuse ISM, its precise composition remains controver- 
sial. Signiflcant clues can be derived from observa- 
tions. The prominent features observed in its spec- 
trum at 9.7 and 18 /im indicate the presence of silicate 
grains. By converse, the origin of the 2175 A feature is 
somewhat debated, though graphite still stands as the 
most attractive solution (see Mathis 1990). The emis- 
sion bands in the 3-13 /im region indicate the presence 
of PAHs (Puget, Leger, & Boulanger 1985), which to- 
gether with small thermally fluctuating grains con- 
tribute also to produce the warm mid-IR cirrus emis- 
sion. Therefore we included these components in our 
diffuse dust model. We adopted the same mixture 
even for MCs, for which the information are scanty, 
only decreasing the fraction of PAH molecules (see 
below). 

The optical properties of silicate and graphite grains 



have been computed by Laor & Draine (1993) for 
spherical shapes using Mie theory, the Rayleigh-Gans 
approximation and geometric optics. In particular we 
used the cross sections of graphite and silicate grains 
computed by B.T. Draine for 81 grain sizes from 0.001 
to 10 fim in logarithmic steps 5 log a = 0.05, and made 
available via anonymous ftp at astro.princeton.edu. 

In order to get the overall radiative properties of 
the dust mixture, the abundance and the size distri- 
bution of grains must be specified. We started from 
the distribution proposed by Draine & Lee (1984, 
henceforth DL), which is tuned on the optical-UV ex- 
tinction law of the galactic diffuse ISM. On the other 
hand DL pointed out that their model displays signif- 
icant discrepancies with extinction data above 2 /im, 
and that these discrepancies are not surprising, since 
observations at longer A tend to sample lines of sight 
through denser clouds, where grain population may 
differ from that present in more diffuse regions. 

A more severe limitation of DL model is that the 
mid-IR (MIR) cirrus emission at wavelengths A <60 
/im is not reproduced (Fig. ||), because grains are 
always large enough to maintain a low temperature 
thermal equilibrium in the general ISRF. This emis- 
sion requires reprocessing of the field by particles, 
very small grains and/or PAH molecules, reaching 
temperatures higher than those attained with equi- 
librium heating (e.g. Puget et al. 1985; Draine & An- 
derson 1985; Dwek et al. 1997). This requirement 
can be matched by decreasing the lower limit of the 
graphite grain distribution. Since these grains must 
be more numerous than the simple extrapolation of 
the DL power law, a break to a steeper power law is 
introduced below a;,. These adjustments to the DL 
model enhance the MIR cirrus emission, but tend to 
degrade the agreement with the observed extinction 
law in the optical-UV region. However a reasonable 
compromise can be found. For graphite grains we 
adopted the following size distribution: 

dn., \ A^hh a'^i ii at < a < a,„ax , 

da \ Ai uh a^^^'^^a^^ if Omin 

with amin = 8 A, flmax = 0.25 /Xm, flb — 50 A, 

(3i = -3.5, 132 = -4.0 and Ag = 10-25.22 cm^ VH. As 
for silicate grains, we maintained the same size distri- 
bution as DL. As a result, 282 atoms per million H 
(ppM hereafter) of C are locked in dust grains. 

In the code the size distributions have been dis- 
cretized in 20 logarithmic bins following the prescrip- 
tions given by Draine & Malhotra (1993). Once the 
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bathing radiation field is specified, the emissivity of 
grains with radius a > 100 A is computed assum- 
ing that they achieve thermal equilibrium, so that all 
grains of a given composition and radius emit as gray 
bodies at a single temperature. Below this limit a 
temperature distribution for each size bin and compo- 
sition is computed following Guhathakurta & Draine 
(1989), and then the emissivity is obtained by inte- 
grating over this distribution. 

To model the observed IR emission occurring in 
bands at 3.3, 6.2, 7.7, 8.6 and 11.3 /im we include 
also a population of planar PAHs. The optical-UV 
absorption cross-section upah is reported in Fig. |^. 
This has been obtained by averaging the cross sec- 
tions of 6 different PAH mixtures measured by Leger 
et al. (1989) above 1200 A, and smoothly joining the 
resulting mean curve to the coronene cross section 
below this limit. The same cross section is used to 
take into account the effects of PAH on the extinction 
curve, assuming a vanishing albedo (Fig. |^). Desert, 
Boulanger, & Puget (1990) have proposed an analyti- 
cal description of <jpah derived from the observed in- 
terstellar extinction curve in the 1200-3300 A range, 
under the assumption that these particles dominate 
its EUV rise. On the other hand the extreme UV 
(EUV) rise can be ascribed to small grains. Actu- 
ally using DL optical efficiencies, the EUV part of 
the extinction curve is produced by small grains. Our 
(ypAH rests on laboratory measurements and predicts, 
in the local interstellar radiation field, a PAH absorp- 
tion and subsequent IR emission lower by a factor 
^ 0.5 with respect to the Desert et al. (1990) analyt- 
ical approximation. 

In order to compute the IR emission of PAHs, their 
heat capacity Cpah{T) must be specified. A numer- 
ically convenient and accurate enough representation 
of the estimate given by Leger, d'Hendecourt, & De 
Fourneau (1989) is: 

^ ,rj.-. [ 9.25x1G-4t ffr< 800 K, 

' = <^ 2 X 10-4 T + .58 ff 800 < T < 2100 K , 

Craax ^ ^ T > 2100 K , 

(7) 

where C^ax = 3[iVt - 2]fc and Nt = Nc + Nh is 
the total number of atoms (carbon and hydrogen) in 
the molecule. We adopt a population of PAHs with 
a distribution dn/dNc oc iV^^'^^ from Nc = 20 to 
Nc = 280. Smaller molecules are easily destroyed by 
UV photons (Omont 1986). This distribution is quite 
similar to those adopted by other authors (Dwek et 



al. 1997, Desert et al. 1990). Astrophysical PAHs are 
thought to be partially dehydrogenated: probably due 
to the large UV fiux in the emission regions, some of 
the CH bonds, responsible for the 3.3, 8.6 and 11.3 
fim bands, are broken. The number of H atoms in 
a molecule is thus written as Nh = XhNs^h, where 
Ns.H is the number of hydrogen sites and xh is the 
H coverage. The relationship between Ns^h and Nc 
depends on the arrangements of hexagonal cycles in 
the molecule. The ratio Nc/Ns,h is rather constant 
around 2 for laboratory PAHs molecules, for which 
Nc ^ 50. However, the structure of typical inter- 
stellar PAHs is likely closer to that of catacondensed 
PAHs (Omont 1986), which are the most compact and 
stable ones and have the general formula Cgp2_ffgp, i.e. 
Ns,H — (6 7Vc)°'^- For these quasi-circular molecules, 
usually assumed to represent the gross features of in- 
terstellar PAHs, the radius is given by a = 0.9 ^/Nc 
A if Nc >> 1. With the above relationship between 
Ns.H and Nc the typical fiux ratios observed in the 
ISM between CH and CC bands (Mattila et al. 1996) 
is fairly well reproduced by our PAH model setting 
XH = 0.2 

PAHs emissivity was then computed following sub- 
stantially the guidelines given by Xu & De Zotti 
(1989)0. The adopted abundances of PAH molecules 
in the diffuse ISM and in the MCs implies that 18 
and 1.8 ppM of C are locked in this component re- 
spectively. Indeed there are indications showing that 
in denser environments and/or in stronger UV radi- 
ation field the relative number of small particles is 
significantly diminished (Puget & Leger 1989; Kim 
& Martin 1996 and references therein). In particular 
Xu & De Zotti (1989) concluded that PAH molecules 
are less abundant by a factor 10 in our galaxy star 
forming regions than in the diffuse gas. 

As apparent from Fig. |^, this model reproduces 
reasonably well both the extinction from IR to UV 
and the whole cirrus emission. The model predicts 
an extinction below a few observational estimates 
at A ^ 300/im by a factor ~ 10. To account for 

this, Rowan-Robinson (1986, 1992) introduced ad 
hoc modifications of the grain optical properties in 
his discretizcd models. We avoid in general this ap- 

^ Apart two imprecisions in their equations: the factor 2 before 
(Tpjxjj in their equation (13) is wrong since their cross section 
already takes into account the two molecule surfaces. On the 
other hand the integration over solid angle yields a factor 47r 
before in the same equation, which in conclusion should be 
multiplied by a factor 2it on the rhs 
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proach (but see Sec. |3.lD mainly because these data 
have large uncertainties (DL), as it is apparent from 
their scatter. Moreover they refer to dusty environ- 
ments for many respects different from the diffuse 
dust, which is instead the dominant contributor to 
the emission at A 100 fim in our galaxy models. 
Also it has been suggested that the silicate grain ab- 
sorption in the sub-mm wavelength range may decline 
less steeply than (Agladze et al. 1996). 

The adopted mixture, even taking into account the 
contribution of PAH emission features, tends to un- 
derpredict the MIR cirrus emission, in particular the 
shorter A data. A larger quantity of small grains 
and/or PAH molecules would improve the match with 
DIRBE observations, but would also produce a large 
disagreement with observed UV-extinction (Dwek et 
al. 1997). We prefer to adopt a compromise more bal- 
anced toward the extinction law, since in our models 
and in most interesting cases the MIR emission is in 
general dominated by warm dust in MCs rather than 
thermally fluctuating grains in the cirrus, whilst a 
good specification of the extinction properties is cru- 
cial. Moreover, DIRBE cirrus data at short wave- 
length could be affected by significant systematic ef- 
fects (Dwek et al. 1997). 

The total (i.e. cirrus -I- MCs) dust content of a 
galaxy ISM depends on the residual gas and on dust- 
to-gas ratio S. The residual gas mass is provided at 
each time step by the chemical evolution model (see 
above). Determinations of S in our Galaxy ISM and 
in nearby objects range typically from 1/100 to 1/400. 
For the purpose of this paper, where the model is com- 
pared only to galaxies in the local universe, we simply 
set d — 1/110, the standard value of DL model. How- 
ever the code is thought to interpret data on objects 
in very different evolutionary stages, thus we need in 
general a recipe to scale S when the chemical condi- 
tions in the ISM where very different. We adopt the 
simplest assumption d (x Z, with the proportionality 
constant adjusted to have S = 1/110 for Z — Zq. 

A recent difficulty for models relying on popula- 
tions of carbon-based grains, the so-called carbon 
crisis, is connected with the quite uncertain determi- 
nation of cosmic abundance of this element. Indeed 
these models were developed under the hypothesis 
that the measured solar system abundances (either 
in meteorites or solar photosphere) were also typical 
for the ISM. For instance, DL as well as our model, 
requires 282 ppM of C in grains, assuming a graphite 
density of 2.26 gr cm^''. This is 55 % of the solar 



abundance which was adopted by DL as cosmic ref- 
erence abundance (but ~ 80 % of the more recent 
estimate in the solar system ~ 365 ppM quoted by 
Anders & Grevesse 1989). However in the last few 
years evidence has been accumulated that solar abun- 
dances may be not representative of the ISM. In par- 
ticular in young star photospheres carbon is found to 
be only ~ 200 - 250 ppM (e.g. Snow & Witt 1995, 
1996; Gies & Lambert 1992). This, if really represen- 
tative of the ISM abundance, would put restrictions 
difficult to meet by carbon-based dust models (Sofia, 
CardeUi, & Savage 1994; CardeUi et al. 1996; Mathis 
1996), especially when taking into account that about 
140 ppM seem to be in gas phase (CardeUi et al. 1996), 
leaving for dust only about 1 /3 of the C required by 
the DL model. 

Possible ways out have been suggested: for in- 
stance big grains could be porous, likely because they 
are built up by sticking of smaller ones (Mathis & 
Whiffen 1989; Mathis 1996), or dust could be to some 
extent prevented to be incorporated into stars during 
their formation (Snow & Witt 1996). Also the intro- 
duction of carbonaceous grains in a form other than 
graphitic (amorphous carbon, possibly hydrogenated, 
or composite grains), the use of more complex size 
distribution than a simple power law, and an elon- 
gated grain shape tend to alleviate the C requirements 
(Mathis 1996; Kim & Martin, 1996). 

Since these problems are still far from being clar- 
ified and our major concern is to obtain a good de- 
scription of the absorption and emission properties of 
dust, we maintain by now a DL-type model, look- 
ing only for the minimum modifications required to 
obtain an acceptable overall match of the extinction 
curve from the UV to the far-IR and of the galactic 
cirrus emission. It is worth noticing that in the future 
this goal should be achieved with a more efficient use 
of the available cosmic carbon. The main effect of 
this on our spectrophotometric galaxy models will be 
a possible decrease of the required dust-to-gas ratio 
by a factor up to ^ 2. 

2.5. The radiative transfer in dusty media 

As already remarked, to get the SED from 0.1 to 
1000 /zm, we have to solve the transfer equation for 
the radiation in presence of dust in the case of molec- 
ular clouds and diffuse ISM. 
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2.5.1. Radiation transfer in MCs and emerging spec- 
tra 

In the Galaxy virtually all the star formation ac- 
tivity resides in molecular clouds. Maps of CO and 
other tracers, show that MCs are non-uniform, highly 
structured objects, containing density-enhanced re- 
gions, the cores, wherein star formation actually oc- 
curs. This implies the clustering of young stars at 
different locations inside giant molecular clouds, as 
confirmed by infrared imaging. The first evolution- 
ary stages are hidden at optical wavelengths and a 
significant fraction, if not all, of the energy of young 
stellar objects is reprocessed by dust and radiated in 
the IR, only a minor fraction being reradiated as re- 
combination lines. The powerful stellar winds and 
outflows, and the ionizing flux from massive stars, all 
contribute to the destruction of the molecular clouds 
in a time scale comparable to the lifetime of OB stars, 
~ 10^ — 10^ yr. Thus stars gradually get rid of their 
parent gas and become visible at optical wavelengths. 
The typical dust densities in these objects are so high 
that even IR photons are self absorbed, thus to com- 
pute the emitted SED the radiative transfer equation 
must be solved. 

The complex evolution depicted above is simulated 
as follows. A fraction /,„c of the model gas mass Mgas 
at tQ is ascribed to the dense phase under discussion 
in this section. Recent estimates in our Galaxy sug- 
gest that half of the hydrogen mass is molecular H2, 
mainly in clouds with diameters ^ 10 pc. The molec- 
ular gas Mmc — fmc Afgas IS then sub-divided into 
spherical clouds of assigned mass and radius, nimc 
and Tmc, which may range in the typical observed in- 
tervals ~ 10^ — 10^ Mq and ~ 10 — 50 pc respectively. 
It is then supposed that each generation of stars, rep- 
resented in our scheme by a SSP, is born within the 
cloud and progressively escapes it. This is mimicked 
by linearly decreasing the fraction / of SSP energy 
radiated inside the cloud with its age t: 

r 1 iit<to, 
f= I 2 - t/to if < t < 2 , (8) 
[0 if <> 2 . 

The model parameter to introduced by this relation 
sets the fraction of starlight that can escape the star- 
bursting region. The starlight locked up inside the 
cloud is approximated as a single central source. The 
cloud optical depth is fixed by the cloud mass and 
radius and by the dust-to-gas ratio. The emerging 



spectrum is obtained by solving the radiative transfer 
through the cloud with the code described by Granato 
& Danese (1994). 

The distribution of stars in real GMCs implies a 
dust temperature distribution with many hot spots 
and cooler regions randomly distributed. A complete 
discussion of the effects of different approximations 
in this complex geometrical situation, including the 
single central source adopted here, can be found in 
Kriigel & Siebenmorgen (1994, see in particular their 
Fig. la). An obvious drawback of our approach is an 
overestimate of the amount of very hot dust around 
the source, with respect to the more realistic approx- 
imation in which the stars are split in many sources 
at different locations. On the other hand the approx- 
imation of a central point source is not much different 
from the IR spectrum predicted by a full treatment 
with hot spots, if the emission from the whole cloud 
is considered (Kriigel & Siebenmorgen 1994). 

Since the treatment of many hot spots in the cloud 
would introduce many other geometrical parameters 
and would slow down considerably our code due to 
the loss of symmetry, we simply treat the maximum 
temperature Tg of the grains at the inner edge as a 
parameter, summarizing in some way the geometri- 
cal parameters which would result in a lower 'true' 
maximum temperature: lowering the grain 'sublima- 
tion' temperature produces a lower average tempera- 
ture with a less opaque cloud. This brings the overall 
spectrum very close to that predicted with a proper 
treatment of the hot spots. It is anyway interesting 
to note that Gordon et al. (1997) found that their 
shell geometry, similar to that adopted here, is suited 
to explain the observed optical and UV properties of 
star burst galaxies, while the dusty geometry, in which 
dust and stars share the same spatial distribution, is 
not, even allowing for dumpiness. 

For typical values of the relevant parameters our 
model predicts IR spectra of star forming regions 
peaking around 40-60 /im (in a vL^ plot), depend- 
ing on the optical depth, with a rather steep decrease 
at longer wavelengths. We found that this behavior 
reproduces quite well the typical range of IR spectra 
of galactic YSOs, H II, and star forming regions (e.g. 
Rowan-Robinson 1979; Ward-Thompson & Robson 
1990; Men'shchikov & Henning 1996). For instance 
in Fig. I the observed SED of W49A, a huge star- 
forming region composed of H II regions and giant 
molecular clouds, is compared with the spectrum of 
the MC component we use in the fit of the overall 
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SED of Arp 220. Indeed we think that this ultralu- 
minous infrared galaxy is a good benchmark for our 
MCs model, having an IR spectrum dominated by 
this component (see Section . 



2.5.2. Propagation in the diffuse ISM and cirrus 
emission 

Before escaping the galaxy the light arising from 
stellar populations and from molecular clouds inter- 
acts with the diffuse dust component. 

We adopt a simplified treatment of radiative trans- 
fer in the diffuse gas which ignores dust self-absorption 
and approximates the effects of optical-UV scattering 
by means of an effective optical depth, given by the 
geometrical mean of the absorption and scattering ef- 
ficiencies Tgyj = Ta{Ta + Tg) (Rybicky & Lightman 
1979, p. 36). Indeed the relatively low opacity of 
dust in the IR regime, where the dust emission oc- 
curs, implies that in most cases of interest the dif- 
fuse ISM is transparent to its own photons. Our ap- 
proximation for combined scattering and absorption 
processes is rigorously applicable only to an infinite 
homogeneous medium and isotropic scattering. How- 
ever we checked, by comparing our results with those 
obtained by Witt et al. (1992) by means of a Mon- 
tecarlo radiative transfer code including anisotropic 
scattering, that it is fairly good in most "real-world" 
geometrical arrangements (see Fig. ^ 

The galaxy is subdivided into small volume ele- 
ments Vi. The local (angle averaged) radiation field 
in the i-th element due to the extinguished emissions 
of free stars and MCs from all the elements is com- 
puted from: 



Vk [j 



mc 
\,k 



k 



(9) 



where Teff^\{i,k) and r{i,k) are the effective optical 
thickness and the distance between the elements i and 
k respectively. Then the local dus t emissivity is calcu- 
lated as described in section 2.4, Finally the specific 



flux measured by an external observer in a given di- 
rection 6 is derived as a sum over the galaxy of the 
extinguished emissivity of free stars, MCs, and diffuse 
dust: 

Fx{e) - A^Y. ^kix,k CXp(-Te//,A(fc,^)) , (10) 



where Teff.\{k,9) is the optical thickness from the 
element k to the outskirts of the galaxy along the 
direction 9. 



In the code the spheroidal systems are assumed to 
extend up to the tidal radius rt — 10^'^rc (of the star 
component), while exponential disks are truncated at 
6Rd, where Rd is the largest of the three components. 
Since more than 98% of the mass is included within 
these radii, adoption of larger cut off radii would not 
produce any significant difference. Deviations from 
the energy balance between dust emission and absorp- 
tion, which mainly depend on the number of volume 
elements in which the galaxy is subdivided, are kept 
within < 2 - 5%. 

3. Comparison with the data 

A few examples, showing the capability of our 
model to reproduce the data on starbursting and nor- 
mal galaxies in the local universe, are discussed below. 
We restrict ourselves to galaxies for which available 
photometric observations allow a precise evaluation 
of the SED for UV to far-IR. The fitting parameters 
are reported in Tables |l| and ^ for the chemical model 
(step 1) and for the photometric model (step 2) re- 
spectively. Table ^ summarizes a few relevant quanti- 
ties derived from those parameters. Since the purpose 
of this paper is to present a procedure to take into ac- 
count the radiative effects of the ISM, the parameters 
tinf and v (Table |l|), which rule the SF history of the 
galaxy, where set following the general results of pa- 
pers dealing with the chemical evolution of different 
types of galaxies (see Matteucci 1996 for a review), 
under the major constraint to get a suitable amount 
of residual gas. The following parameters have been 
fixed to reasonable values: galactic age to = 13 Gyr; 
dust-to-gas mass ratio S — 9 x 10~^; exponent of 
Schimdt SFR fc = 1; maximum temperature of dust in 
MCs Ts = 400 K, mass of single MC rumc = IO^Mq. 
Moreover the mass limits of the Salpeter IMF are the 
standard M„p = lOOM© and Miow = 0.1, but for 
starburst galaxies where the estimates of dynamical 



mass require Mu 



0.2, as discussed below. We also 



remind that, in order to keep the number of free pa- 
rameters to the minimum required by present quality 
data, we set the scale lengths of stars and gas distri- 
butions in spirals to the same values, i.e. i?^ = 
and z*^z^. 

As previously noticed, the SED of the MCs com- 
ponent is set mainly by its optical depth (Table ||) 
oc nimc/rmci which is in some sense the true "fit- 
ting parameter" for this component (together with 
to), rather than r„jc reported in Table 0. In other 
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words the values found for r^c depend on having set 
rrimc = W^Mq, as typical for Giant Molecular Clouds 
in the Galaxy, and r^c could well be varied by a factor 
~ 2 producing almost equivalent fits, provided nimc 
would be adjusted to keep nimc/r'^c constant. 

3.1. Starburst galaixies 

M82 

In the prototype starburst galaxy M82 the burst 
was probably triggered by the interaction with M81, 
some 10^ yr ago (Solinger, Morrison, & Markert 
1977). Thanks to the proximity of this system {D = 
3.25 Mpc) a wealth of data do exist, providing a well 
sampled full coverage of the SED at different angu- 
lar resolutions, as well as other observational con- 
straints. The fit (Fig. ^ is obtained by evolving for 
13 Gyr an open system with a final baryonic mass 
of 1.8 X 10^° M0, a factor 6 less than the estimated 
total dynamical mass (Doane & Mathews 1993). We 
adopted Mi = 0.2 Mq as lower limit of the Salpeter 
IMF. A lower value would require a baryonic mass 
closer to the dynamical mass. The assumed parame- 
ters provide a SFR raising from to about 3 Mq /yr 
in the first 3 Gyr, then smoothly declining to 1.35 
Mo/yr at t = 13 Gyr. To this gentle star formation 
history, which leaves a gas fraction of 0.064, we have 
superposed an exponential burst processing 18% of 
the residual gas in the last 5x lO^r, with an e-folding 
time also of 5 x 10^ yr (Figure 0). Note that our re- 
sults are almost independent on the precise evolution 
of the SFR in the burst: even a constant burst yields 
very similar estimates of the relevant quantities. Also 
the duration of the burst could be doubled or halved 
with small adjustment of the other parameters with- 
out damaging the quality of the SED fit. Thus we end 
with a total gas mass of 8.6 x 10* Mq, 8% of which 
is ascribed to the molecular component, organized in 
clouds with TOmc = 10^ Mq and rmc — 16 pc. These 
clouds reprocess almost completely the starlight due 
to the burst. The system is assumed to follow a King 
profile with r* = 150 pc for stars and = 200 for 
diffuse gas. 

The masses we ascribe to the various components 
favorably compare with radio estimates of total gas 
masses - IO^Mq (Solinger et al. 1977), with CO 
determinations of gas in molecular form ~ 10* Mq 
(Lo et al. 1987; Wild et al. 1992), and with an up- 
per limit ^ 3 x 10* Mq to the mass of stars formed 

in the burst, suggested by dynamical considerations 



(McLeod et al. 1993). Also, in our proposed model 
the predicted supernovae rate (SNR) is between 0.05 
and 0.1 yr~^, depending on the adopted lower limit 
of stellar mass yielding to supernovae explosion 9 or 6 
A/0 respectively. Observational estimates suggest a 
SNR in the range 0.07-0.3 yr"! (McLeod et al. 1993; 
Doane & Mathews 1993). Increasing the lower limit 
of the IMF in the burst to, for instance. Mi = 1.0 
would increase the SNR by a factor ~ 2. In conclusion 
the model required to nicely reproduce the observed 
SED turns out to be in agreement with independent 
estimates of masses and of SNR. 

NGC 6090 

This strongly interacting galaxy at 175 Mpc (Hq = 
50 km/s/Mpc) has been observed by ISO from 2.5 
to 200 tJLui (Acosta-Pulido et al. 1996). The SED 
resulting from the combination of these data with 
previously published optical photometry is nicely re- 
produced by our model (Fig. ^). The differences in 
the parameters with respect to M82, apart from an 
up-scale in involved baryonic mass (4.1 x IQ^^ Mq), 
are aimed at enhancing the diffuse dust emission at 
A ^ 100 /im, which, as noticed by Acosta-Pulido et 
al., is not reproduced by pubhshed starburst models. 
Thus the pre-burst SFR has been adjusted to leave a 
larger gas fraction (0.127), 3.2 % of which is processed 
by the burst. The molecular star- forming clouds ac- 
counts only for 0.5 % of the gas left by the burst. 

Arp 220 

Arp 220 is an archetypal ultraluminous infrared 
galaxy (ULIRG), most likely the result of a recent 
merging between two gas-rich galaxies, for which 
ISOPHOT data have been now published (Klaas et al. 
1997), and whose K-band light profile resembles that 
of a typical elliptical galaxy. In this object there is 
evidence of both starburst as well as Seyfert activity, 
but ISO spectroscopy led to the conclusion that the 
IR luminosity is primarily ( ^ 90 %) powered by star- 
burst (Sturm et al. 1996). The fitting model requires 
a strong burst, converting into stars 2.5 x 10^° M©, 
i.e. as much as 11% of the total baryonic mass, and 
a high fraction 50% of residual gas in star-forming 
clouds. The gas fraction before the burst is 0.25. As 
a result, the IR and sub-millimetric emission is every- 
where dominated by dust associated with star forming 
regions. The mass in molecular gas is ^ 1.6 x 10^*^ Mq, 
in agreement with CO estimates (Solomon et al. 1997; 
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Scoville, Yun, & Bryant 1997). From a determination 
of the CO rotation curve Scoville et al. (1997) infer a 
dynamical mass of 3 — 6 x IO^^A/q within r ~ 1.5 
kpc (depending on the adopted potential), in agree- 
ment with the mass requested by our model within 
the same radius ~ 4.1 x lO^^M©. 

We find that the sub-millimetric data of Arp 220 
are underpredicted by a factor ~ 2 — 3 by any model 
which successfully reproduces the remaining SED, un- 
less the long wavelength decline of grain absorption 
efficiency is reduced, say from cx to cx A"""^'^ 
above ~ lOO/itm (see Agladze et al. 1996 for laboratory 
measurements supporting this possibility for silicate 
grains). In Fig. |^ we report fits obtained adopting 
either the modified as well as the standard DL cross- 
sections. 

3.2. Normal galaxies 

In the following we examine objects that do not 
show clear signs of enhanced SF or nuclear activities 
in the SED nor in the morphology, namely 3 late- 
type spirals and a giant elliptical template. Indeed 
elliptical galaxies form a relatively homogeneous class, 
while this is not the case for spirals. We fit spirals 
later than Sb, which are clearly disk-dominated, and 
we adopt for them an exponential geometry (see Sec. 

II)- 

M51 

In Fig. |l^ we present a fit to the SED of the nearly 
face-on {i = 20°, TuUy 1974) Sbc galaxy M51 (NGC 
5194), taken at a distance D = 9.6 Mpc (Sandage & 
Tammann 1975). Since the discrepancy between ISO 
and IRAS data at 60 and 100 /im could be due to an 
excess of the ISO point spread function (Hippelein et 
al. 1996), we fit the IRAS data. The model has a bary- 
onic mass of 1.55 x 10^^ Mq, whose evolution leaves at 
13 Gyr a gas fraction of 0.067, 70% of which in molec- 
ular form. The mass of ISM in the molecular and 
diffuse components, as well as their relative fraction 
is in agreement, within a factor < 2, with estimates 
by Scoville & Young (1983), Young et al. (1989), and 
Devereux & Young (1990). Gas and stars are ex- 
ponentially distributed with the same scale-lengths, 
4.7 and 0.4 kpc for the radial and vertical scales re- 
spectively. The former value is in good agreement 
with those estimated by Beckman et al. (1996) from 
observed brightness distributions in different optical 
bands. Due to the low energy output observed at 



A ;S 0.2/im, stars born during the last 5 x 10^ yr are 
hidden inside molecular clouds. 

The striking correspondence between the 15 ^m 
emission mapped by ISOCAM and the emission 
indicates that the MIR is powered by recent star- 
formation (Sauvage et al. 1996). As for the far-IR 
(FIR), somewhat contradictory claims have been re- 
ported: Devereux & Young (1992), comparing the ra- 
dial distributions of FIR, Hq,, H I, and II2 emission, 
concluded that the same holds true in the range 40- 
1000 fim, while Hippelein et al. (1996) found no obvi- 
ous correlation between FIR ISOPHOT maps and 
fluxes. This apparently complex situation is not sur- 
prising, since according to the model the MIR emis- 
sion is provided by the MCs component, while above 
~ 60 /im the diffuse dust gives a comparable contri- 
bution, which however is in part powered by young 
stars. 

Ml 00 

The Sbc galaxy MlOO (NGC 4321) is the largest 
spiral in the Virgo Cluster. We adopt a distance 
D = 20 Mpc and an inclination angle i — 30°. The 
model in Fig. |ll| has a baryonic mass of 2 x 10^^ Mq, 
with a residual gas fraction of 0.048, 80% in molecular 
form. Indeed, according to the gas masses estimated 
by Young et al. (1989) and Devereux & Young (1990), 
larger than ours by a factor ~ 2, the ISM seems to 
be dominated by the molecular component. The lo- 
cation of MlOO in the central regions of the Virgo 
Cluster could affect some of its properties, as the H 
I distribution that shows a sharp edge in correspon- 
dence to the optical radius (Knapcn et al. 1993). This 
could be the reason for the warm FIR SED observed 
in MlOO, as compared to M51 and NGC 6946, de- 
spite their morphological similarity. The steep decline 
of the spectrum in the sub-mm region has been in- 
terpreted by Stark et al. (1989) as indicating that the 
emitting grains are warm and small. In our model this 
is not required since the fit is obtained with standard 
diffuse dust. Stars and dust share the same radial 
and vertical scale-lengths, 5 and 0.4 kpc respectively, 
in agreement with the values derived by Beckman et 
al. (1996) for the stellar component and with the H I 
distribution. 

NGC 6946 

The fit to this Scd galaxy at distance D = 6.72 
Mpc (Rice et al. 1988) and inclination angle i = 34° 
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(Considere & Athanassoula 1988), is shown in Fig. [T^ . 
The total gas mass in the model, 1.04 x 10^° Mq i.e. 
8.3% of the baryonic mass, agrees with the molecu- 
lar plus neutral hydrogen mass given by Young et al. 
(1989) or Devereux & Young (1990) (reported at 6.72 
Mpc), while the fraction we ascribe to the molecular 
component is a factor ~ 2 higher. 

The contribution of young and old stellar popula- 
tions to dust heating in this object has been consid- 
ered by several authors. Devereux & Young (1993), 
comparing the radial distributions of FIR, Hq,, H I, 
and H2, conclude that the 40 to 1000 /xm luminosity 
is dominated by dust associated with molecular gas 
heated by young stars. However ISO data reported 
by Tuffs ct al. (1996) reveal an extended FIR emission 
out to a radius of 8', with a scale-length similar to 
the R-band one, while little Hq, emission has been de- 
tected beyond r ~ 6'. This suggests that part of the 
observed cold FIR SED of NGC 6946 could be due to 
dust associated with a diffuse H I gas, which indeed 
extends out to r '--^ 15' (Boulanger & Viallefond 1992), 
and heated mainly by old stellar populations. Malho- 
tra et al. (1996) find that the radial scale-lengths at 7 
and 15 /im are similar to those in Hq, and H2 but much 
shorter than those in R-band and H I, consistent with 
a warm dust emission primarily heated by young mas- 
sive stars. In the proposed fit the MIR emission arises 
from MCs dust heated by newly born stars. The SED 
above ^ 60 /im is almost equally contributed by MCs 
and diffuse dust, but, since to — 2.5 Myr, even the 
latter is predominantly heated by young stars. 

Tacconi & Young (1986) estimate a scale-length of 
9.6 kpc for the cold ISM component, whilst for the 
star distribution their values range from 4 to 8 kpc, 
depending on the band and on the galactic component 
(disc or arms) considered. The model scale-lengths 
are set to 8 for both stars and the diffuse gas. The 
vertical scale-lengths are 1 kpc to keep the diffuse 
dust emission sufficiently cold. 

Giant Ellipticals 

Ellipticals constitute a class of objects with fairly 
homogeneous spectral properties. It is then interest- 
ing to test our model against an average SED of gi- 
ant E galaxies, constructed combining the Arimoto's 
template (1996) from 0.12 to 2.2 /xm with the me- 
dian IRAS over B band fluxes of bright ellipticals es- 
timated by Mazzei et al. (1994) and the average (K-L) 
color computed by Impey et al. (1986) (Fig. |l^). The 
fit has been obtained with a 'classical' star forma- 



tion history model for ellipticals: an open model with 
infall and galactic wind. The high efficiency v = 2 
causes a huge SFR in the first 1.2 Gyr interrupted 
by the onset of the galactic wind (see Tantalo et al. 
1996 for a discussion of these values). The object is 
observed at 13 Gyr, when a comparatively very small 
fraction of galactic mass in a diffuse dusty component 
is sufficient to produce the observed IRAS emission. 
Note that this gas is not provided by the chemical 
model, according to which the bulk of the ISM is 
ejected by the galactic wind which stops the star for- 
mation. Thus in this case the mass in the diffuse dust 
1.5 X 10^ Mq is a fitting parameter of our photometric 
model. This dusty ISM could arise from evolved stars 
of the passively evolving galaxy, by cooling fiows or 
by merging activity. It is also worth noticing that Tc 
needs to be much greater for the ISM than for stars 
(6 and 0.4 kpc respectively): for smaller values of 
the cirrus emission would be too warm because dust 
would be more concentrated in the central regions 
where the radiation field is higher. Other authors al- 
ready suggested that diffuse dust in elliptical galaxies 
needs to be less concentrated than stars, on the ba- 
sis of either optical color gradients (e.g. Wise & Silva 
1996) or IRAS colors (e.g. Tsai k Mathews 1996). 
As discussed in section 2.5.2 , an alternative proposed 
way to describe this lower concentration is to adopt 
the same Tc for both stars and diffuse dust, decreas- 
ing the exponent 7 = 1.5 in the King law of stars. 
However for this SED, adjusting 7*^ and Tc with the 
constrain = we found only marginally accept- 
able fits in the IRAS regime, significantly worse than 
that presented in Fig. |l^. 

4. Discussion 

When reproducing starburst galaxies, our model 
is characterized by the variations of 6 parameters re- 
lated to the total mass and star formation history and 
of 5 geometrical parameters (Tables || and For 
normal disc galaxies, the parameters become 3 and 5 
respectively. In order to fit the elliptical galaxy tem- 
plate, 7 parameters have been adjusted, including the 
epoch at which the SF activity is stopped. 

The masses reasonably constrained by matching 
the SED are the baryonic mass (though dependent on 
the IMF, particularly its lower limit) mostly from the 
energy emitted in the near-IR (NIR) by the old stel- 
lar component, and the mass of dust, from its MIR- 
FIR emission. Therefore the mass of gas present in 
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the galaxy is constrained from the fit, provided S is 
known. Our assumption S = 9x 10^^ is likely the best 
guess for our Galaxy, however in different systems and 
in the Galaxy itself, variations of a factor of a few are 
commonly quoted, e.g. Rand, Kulkarni, & Rice (1992) 
for M51 or Tuffs et al. (1996) for NGC 6946. As a con- 
sequence, our estimated gas masses suffer by a similar 
uncertainty. The SFR history of the old stellar com- 
ponent is constrained by the spectrum and, for star- 
burst galaxies, also by the requirement that a certain 
amount of gas is available for the latest burst. The 
three analyzed starbursts exhibit rather different star 
formation histories, with the average SFR of the old 
component inversely proportional to the strength of 
the SFR in the burst. Estimates of masses in stars 
and g well as of SF and SN rates from observa- 
tions other than broad band spectra, discussed in the 
previous sections and below, nicely agree with those 
provided by our SED fitting. 

4.1. Starbursts and ULIRGs 

The SFR in a starbursting region can be estimated 
on the basis of the FIR luminosity, under the assump- 
tion that dust reradiates the overall bolometric lu- 
minosity and after calibration with stellar synthesis 
models. Kennicutt (1997) proposes SFR[Moyr~^] = 
Lpifi/{2.2 X lO'*'^ ergss"^), with calibration derived 
from models of Leitherer & Heckman (1995). This 
relationship properly applies to star forming regions, 
while the extension to a galaxy as a whole is danger- 
ous, since a non negligible portion of the FIR lumi- 
nosity may arise from cirrus emission powered also by 
relatively old stars. However our model accurately de- 
termines the warm component associated to the MCs 
where star formation is occurring (im.c in Table |^). 
For M82 our fit predicts an overall FIR luminosity as- 
sociated to the MCs component of 1.1 x 10*'' ergs/s, 
which following Kennicutt translates in SFR=5.0 Mq 
yr""'. Similarly, the SFRs inferred from the luminos- 
ity of the warm component is 57 Mq yr~^ and 460 
Mq yr-1 for NGC 6090 and Arp 220 respectively. 
These figures are in very good agreement with those 
derived by fitting the SEDs and reported in Table ^. 
Thus while Kennicutt's calibration is confirmed, we 
stress that this refers only to the warm component, 
which contributes ~ 50% ofLpiR in NGC 6090. This 
source of uncertainty adds to the obvious effect of the 
adopted IMF in the conversion from observed IR lu- 
minosity to SFRs. 

Whilst the very recent (say in the last 10 Myr) SFR 



is relatively well constrained by the observed warm IR 
emission, the burst duration, and therefore the total 
mass converted into stars, can be varied within a fac- 
tor ~ 2, still yielding the correct Lpjfj and spectra 
after readjustments of other parameters. However the 
supernovae rates deduced from observations put fur- 
ther constraints on the average SFR and duration of 
the burst. In the well studied case of M82 the model 
predicts a SN rate in the interval deduced by obser- 
vations. 

The SFR can also be inferred from the ioniz- 
ing luminosity L^yc, which can be derived from re- 
combination lines. With the adopted IMF we find 
iiyciW^ergss-i] ~ 0.017SFR(MQj/r-i), which holds 
with very small deviations for the six studied star- 
forming galaxies (see Table The uncertainties 
connected to extinction corrections of observed line 
fluxes are minimized employing transitions occurring 
in the IR regime, such as Br7 line (A = 2.17/im). 
Calzetti (1997) and Kennicutt (1997), using Leitherer 
& Heckmann (1995) results, find SFR[Moyr-i] = 
L{Brj) /1.6 X lO^^ergss^-' adopting a Salpeter IMF 
within 0.1-100 Mq or SFR[Moyr-i] = L{Brj)/3.7x 
lO'^'^ergs s~^ when the mass range is 0.1-30 Mq. 

The (Brj) hnc luminosity L(i3r7) = 9.2 x lO***^ ergs s~^ 
of NGC 6090 (Calzetti, Kinney, & Storchi-Bergmann 
1996) corresponds to a SFR=56 Mq yr~^ with the 
assumption of the 0.1-100 Mq mass range, only 18% 
less than that used in our model, while the smaller 
range in the IMF would predict a SFR higher by 
a factor 4. Kennicutt (1997) has shown that there 
is a clear trend for SFRs derived from Lpjji to be 
larger than those estimated from the L{Brj), suggest- 
ing that extinction is non negligible even at near-IR 
wavelengths. 

Genzel et al. (1997) reported the ratio of the far-IR 
to Lyman continuum luminosity Lpm/ L^yc of star- 
bursts (including M82) and ULIRGs, derived from 
near and mid-IR recombination lines. For 12 star- 
burst galaxies the median value is Lpjn/ L^yc — 16. 
The model ratio depends on the adopted IMF. Our 
fit to M82 overall SED yields Lpm/LLyc ^ 12, where 
L^yc is the unextincted luminosity of the (young) 
stellar populations below 912 A in very good agree- 
ment with the value found by Genzel et al. (1997). 
For NGC 6090 the fit predicts Lpm/LLyc ^ 10, 
using Lmc to evaluate the model Lf/_r powered by 
the burst, well within the range of the values in- 
ferred by the same authors. As for Arp 220, which 
properly belongs to the ULIRG class, higher values 
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15 ;S LfirJ L^yc ^1 112 are inferred from IR recombi- 
nation lines, whereas our model predicts Lpm/ L^yc ~ 
11. The discrepancy is likely due to the uncertain 
large correction for IR extinction in this object. Ac- 
tually Genzel et al. (1997) use a screen obscuration 
with Ay = 45, while we find Ay — 150 in MCs of this 
object. Interestingly enough, the 15 ULIRGs (includ- 
ing Arp 220) studied by Genzel et al. (1997), exhibit 
a median Lpm/ L^yc ~ 40, larger by a factor of ~ 2.5 
than that inferred for starburst galaxies. This result 
can be explained by larger obscurations or by a softer 
intrinsic Lyman continuum, pointing to an IMF less 
rich in massive stars or to an older starburst. In the 
models we used bursts began 0.05 Gyr ago, but still 
rather active due to the large e-folding time te = 0.05 
Gyr. 

It is worth noticing that our model does not ac- 
count for the ionizing radiation converted into recom- 
bination lines, which are much less absorbed by dust 
than the Lyman continuum. In the present version of 
the code the UV radiation can be only converted di- 
rectly to IR photons through dust absorption. How- 
ever the energy budget is not significantly affected, 
since the observed luminosity in recombination lines 
is only a few percent of the bolometric luminosity of 
starburst galaxies (e.g. Genzel et al. 1997). 

The observed UV emission at A > 1000 A IS cru- 
cial in determining the fraction of the SSPs escaped 
from the parent molecular clouds. In our model the 
UV flux emerging from a single galaxy depends on the 
age to after which a SSP starts to get out from the 
molecular cloud. The UV photons are affected also 
by absorption in the diffuse ISM. In the case of M82 
and Arp 220 the measured UV flux is a tiny fraction 
of the bolometric luminosity (see Figg. ^ and This 
implies that the stars of the burst are still inside their 
parent clouds. Indeed in these two objects to is about 
equal to the time since burst ignition 0.05 Gyr. By 
converse for NGC 6090 the data show that the UV 
is not a negligible fraction of the total, implying for 
our model ta — tturst > to = 0.018 Gyr. As a re- 
sult stars with lifetime > 18 Myr (corresponding to 
M <, lOM©) are no more embedded within the clouds 
and contribute significantly to the ISRF. Actually the 
luminosity we ascribe to the diffuse cold component 
equals the warm one in NGC 6090, whereas it is 25% 
and 9% for M82 and Arp 220 respectively. 

Calzetti, Kinney, & Storchi-Bergmann (1994) com- 
puted the "intrinsic depth" of the 2175 A dust absorp- 



tion feature r] (see their Eq. 23) in a number of star- 
burst galaxies, finding that —0.15 < 77 < 0. For the 
three starbursts we find 77 ranging from -0.15 to -0.08. 
Since the dust properties we adopt follows very well 
those of our galaxy, at least in these objects the rel- 
ative weakness of the 2175 A feature can be ascribed 
to the adopted geometry, without invoking possible 
variations of extinction curve. 

In conclusion for starbursts, the spectral coverage 
from UV to submillimeter wavelengths allows a robust 
evaluation of the luminosity of stars involved in the 
bursts and of the general interstellar radiation field. 
The latter is contributed by the long-lived stellar pop- 
ulations and by the stars produced in the burst but 
old enough to get out from the original MCs. Our 
model well describes the complexity of a starburst 
galaxy, and is flexible enough to reproduce the differ- 
ences among them. It uses a reasonable number of 
parameters, which have a well defined physical mean- 
ing and compare favorably to values derived from ob- 
servations other than broad band spectra. 

4.2. Spirals 

The SEDs of three spiral galaxies, namely M51, 
Ml 00, and NGC 6946, have been well reproduced by 
our model. An interesting result is the significantly 
lower to in spirals with respect to starbursting ob- 
jects. This derives from their smaller observed ratios 
Lfir/Luv- In the three spirals analyzed here, the 
newly born SSPs quite soon get out from the parent 
clouds and significantly contribute to the heating of 
the diffuse dust. The expected contribution is larger 
for NGC 6946 and MlOO than for M51. As a conse- 
quence, a significant portion of the far-IR luminosity 
emitted by the diffuse dust (almost all for NGC 6946 
and MlOO) should be inserted in the budget to infer 
the SFR in spiral galaxies, in agreement with the find- 
ings of Devereux et al. (1994) and Buat & Xu (1996). 
Taking into account this, the SFRs derived from the 
FIR luminosities, through the conversion suggested 
by Kennicutt (1997), are in very good agreement with 
those used by the model. 

The SFRs in our models, through the conversion 
of Kennicutt (1997) i(iJa)[ergs s"!] = 1.26 x W^^ x 
SFRIMqjt-^], yield luminosities of 7.6, 8.8, and 
7.610^1 ergs s'^ for M51, MlOO, and NGC 6946 re- 
spectively, to be compared with observed values (Ken- 
nicutt, Tamblyn, & Congdon 1994, reported to our 
adopted distances) of 4.9, 5.0, and 2.2 x 10^^ ergs 
s~^. These lower values would imply an internal 
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extinction of Aho — 0.5, 0.6, and 1.3 magnitudes in 
the three spirals, consistent with the average value 1.1 
adopted by Kennicutt (1997), based on a comparison 
of free-free radio and fluxes. 

It is also worth noticing that masses in stars, dust, 
and gas derived by the model, as well as geometric 
parameters such as the disc scale length are in good 
agreement with independent estimates. Our model 
includes the main physical aspects of star formation, 
stellar evolution, and dust absorption in spiral galax- 
ies and allows a full exploitation of the broad band 
data in order to explain their present status and their 
past history. 

4.3. Ellipticals 

The model with a relatively small number of pa- 
rameters is able to produce a very good fit to the 



template SED of giant elliptical galaxies (see Fig. 13). 
Most excitement about spheroidal galaxies is related 
to their evolution. Indeed it has been suggested that, 
during the initial phases of star formation, they might 
look very similar to local violent starbursts (such as 
Arp 220), since the chemical enrichment and, as a 
consequence, the dust formation are very quick pro- 
cesses when SFRs are very high (Mazzei et al. 1994; 
Franceschini et al. 1994). This initial phase, if con- 
fined to high enough redshifts z < 2 — 2.5, is the 
most natural way to produce the total energy and 
the shape of the Far-IR Background (Franceschini 
et al. 1991; Franceschini et al. 1994; Burigana et al. 
1997), which has been tentatively detected by Puget 
et al. (1996) and recently confirmed by Hauser et al. 
(1998) and Fixsen et al. (1998). Also we expect that 
in the ISO extra deep surveys we may detect individ- 
ual dusty spheroidal galaxies at substantial redshifts. 
Our model is able to describe in detail the evolution 
of elliptical galaxies in all the relevant range of fre- 
quencies, and will be exploited to a comprehensive 
discussion of galaxy counts, related statistics, and as- 
trophysical backgrounds. 



ties such as SFR, the IMF and the past history of the 
galaxies. The model can be implemented in studies 
of local starburst and normal galaxies. It can also be 
used to trace back the history of the different classes 
of galaxies and, by confronting predictions and obser- 
vations of number counts and redshift distributions in 
different bands, to understand their cosmological evo- 
lution. In particular we will use the model to investi- 
gate the role of dust during the early phase of galaxy 
formation, exploiting available and soon coming data 
on galaxy counts at wavelengths ranging from UV to 
radio bands. 

Updated information on the code described in 
this study (GRASIL) can be found on th e WEB at 
http:/ / asterix.pd. astro, it/homepage/ gian/ . 
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4.4. Conclusion 

We showed that the model and the related numer- 
ical code are extremely efficient in deriving a wealth 
of information from broad band spectra of starburst 
and spiral galaxies. Indeed, when SEDs from the UV 
to the sub-mm range are available, we are able to 
quantify the effects of dust reprocessing on observa- 
tions, gaining information on very substantial quanti- 
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Fig. 1. — Scheme of the components included in our 
model computations and their adopted geometry. 
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This 2-column preprint was prepared with the AAS lATJiiX 
macros v4.0. 



Fig. 2. — Solid line: cross-section of PAH per carbon 
atom, derived from laboratory measurements taken 
from Leger et al. (1989); dashed line: the analytical 
law assumed by Desert et al. (1990) for molecule with 
Nc = 50 C-atoms {Nc controls the cut-off above 0.3 
/im). 
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Fig. 3. — Upper panel: the extinction curves of the 
dust model we adopt for diffuse and molecular gas, 
and that by Draine & Lee (1984) are compared with 
available data. Lower panel: the predicted emissivity 
of grains (ergs cm~^ st~^ ^h^) in the local in- 
terstellar radiation field (Mathis, Mezger, & Panagia 
1983) compared with observations toward the galac- 
tic pole. The orizontal line around 12/im marks the 
flux level obtained by convolving our expected SED 
with the 12/im IRAS passband. Triangles in the 
lower panel are data from Dwek et al. (1997), for 
references to other observations see Rowan-Robinson 
(1992, 1986) 
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Fig. 4.— The observed SED of W49A (Ward- 
Thompson & Robson 1990), a huge star-forming re- 
gion composed of H II regions and GMCs, is compared 
with the spectrum of the MC component we use in the 
fit of the overall SED of Arp 220. 
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Fig. 5. — F/Fo is the ratio between the dust ex- 
tinguished flux and that expected if dust were not 
present. The points have been computed by Witt et 
ah (1992) with a Montecarlo radiative transfer code 
including anisotropic scattering, the dashed hne rep- 
resents the result of our code taking into account only 
absorption, whilst the solid line includes scattering 
with the effective optical depth r^^j = Ta{Ta + Ts). 
The adopted geometry is that defined by Witt et al. 
as elliptical galaxy. 
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Fig. 6.— Fit to the SED of M82. Data are from 
Code k Welch (1982), Soifer et al. (1987), Klein et al. 
(1988), Cohen & Volk (1989), Van Driel et al. (1993), 
Ichicawaet al. (1994,1995). 
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Fig. 7. — M82: SFR as a function of galactic time. 
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Fig. 9. — Arp 220: in this case the wavelength 
dependence of grain cross-section has been modi- 
fied above 100/xm from cx to cx A~^'®. The 
dashed line above 100 /im represents the model pre- 
diction with standard decline. Data are from 
Carico et al. (1988), Sanders et al. (1988), Smith et 
al. (1989), Carico et al. (1990), Wynn-WiUiams & 
Becklin (1993), Rigopoulou et al. (1996), Klaas et al. 
(1997). 



Fig. 8.— Fit to the SED of NGC 6090. Data are 
from Mazzarella & Boroson (1993), Acosta-Pulido et 
al. (1996), Gordon et al. (1997). 
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Table 1 
sf history parameters 



object 


D 


Mg 


1/ 




Mhurst 


^burst 


te 




Mpc 




[Gyr-i] 


[Gyr] 


[IQI^Mq] 


[Gyr] 


[Gyr] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


M82 


3.25 


1.8 


1.2 


9 


0.02 


12.95 


0.05 


NGC 6090 


175 


41 


0.6 


9 


0.16 


12.95 


0.05 


ARP 220 


115 


23 


0.3 


9 


2.5 


12.95 


0.05 


M51 


9.6 


15.5 


0.6 


4 








MlOO 


20 


20 


0.75 


4 








NGC 6946 


6.7 


12.5 


0.6 


5 








SE 




100 


2.0 


0.1 









Note. — Parameters for the star formation history (Sec. 2.1): (2) adopted distance; (3) final baryonic galaxy 
mass; (4) SF efficiency; (5) infall timescale; (6) gas mass converted into stars during the burst; (7) galaxy age at 
beginning of the burst (when included); (8) e-folding time for SFR in the burst. All models have an age of 13 Gyr. 
In the case of gE the SF has been stopped at 1.15 Gyr. 



Table 2 
Geometric parameters 



object 


fmc 




to 


T 




Rd 


Zd 






[pc] 


[Myr] 


[kpc] 


[kpc] 


[kpc] 


[kpc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


M82 


0.08 


16 


57 


0.15 


0.2 






NGC 6090 


0.005 


17 


18 


0.5 


1.0 






ARP 220 


0.5 


10.6 


50 


0.5 


0.5 






M51 


0.7 


14 


8 






4.7 


0.4 


MlOO 


0.8 


15 


3 






5.0 


0.4 


NGC 6946 


0.6 


14 


2.5 






8.0 


1.0 


gE 








0.4 


6.0 







Note. — Parameters for the photometric model estimated from SED fitting: (2) fraction of residual gas in MCs; (3) radius 
of MC; (4) parameter regulating the escape of young stars from MCs (Eq. (5)-(8) parameters for the spatial distribution 
of stars, MCs, and cirrus (Eqs. |l| and ^). For disks the inclination angles have been taken from the literature (see text) and 
are i = 20° (M51), i = 30° (MlOO), and i = 34° (NGC 6946). 
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Table 3 
Derived quantities 



object 


< SFR > 


Mdust 




Tl 


TB 




Lc 


LLyc 




[Me/yr] 


[10'Mq] 








10** ergs/s 


10'"' ergs/s 


10*'' ergs/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


M82 


5.5 


0.8 


25 


0.62 


1.30 


1.1 


0.25 


.09 


NGC 6090 


68 


45 


24 


0.84 


1.37 


13 


13 


1.3 


ARP 220 


580 


30 


58 


1.67 


2.80 


100 


5.7 


8.7 


M51 


6 


10.4 


33 


0.06 


0.28 


1.01 


0.90 


0.11 


MlOO 


7 


9.6 


29 


0.03 


0.15 


0.84 


0.85 


0.12 


NGC 6946 


6 


10.4 


33 


0.03 


0.12 


0.62 


0.53 


0.10 


gE 




0.15 




0.00 


0.01 




0.02 


0.03 



Note. — A few quantities derived from the models: (2) SFR averaged over the last 5 x 10^ yr, which is the time from burst 
onset for starburst models; (3) total mass in dust in the galaxy; (4) 1 /^m optical thickness of the MCs from the centre; (5) 
1 fim 'average' optical thickness of the model, defined such as (observed flux) = (dust-free fiux x exp(— r)) (6) same as (5) 
but at 0.44 ^m; (7) luminosities of the MCs, (8) of the diffuse dust and (10) in the Lyman continuum before dust absorption. 




Fig. 10.— Fit to the SED of the Sbc galaxy M51. 
Data are from Buat et ah (1989), Evans (1995), de 
Vaucouleurs et al. (1991) (RC3), Code & Welch 
(1982), Young et al. (1989), Rice et al. (1988), Dev- 
ereux & Young (1990,1992), Smith (1982), Hippelein 
et al. (1996) (ISO, triangles). 



Fig. 11.— Fit to the SED of the Sbc galaxy MlOO. 
Data are from Buat et al. (1989), Donas et al. (1987), 
De Jong et al. (1994), Stark et al. (1989), RC3, Dev- 
ereux &Young (1990), Young et al. (1989), Helou et 
al. (1988), Knapp et al. (1987). 
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Fig. 12.— Fit to SED of tlie Scd galaxy NGC 6946. 
Data are from Rifatto et al. (1995), RC3, Engargiola 
(1991), Devereux & Young (1993), Rice et al. (1988), 
Tuffs et al. (1996). 
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Fig. 13. — A fit to a template SED of a giant elliptical 
galaxy. 
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